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Abstract 

We point out some problems with the previously-proposed phase diagram 
of the 1*6 spin models. Consideration of the diagram near to the decoupling 
surface using both exact and approximate arguments suggests a modification 
which remedies these deficiencies. With the aid of a new parametrisation of the 
phase space, we study the models numerically, with results which support our 
conjectures. 
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The discrete Zg spin models describe the behaviour of a collection of spins «Sj = 
exp(i0j), where the 9i are integer multiples of ir/3. These spins live on a square 
two-dimensional lattice and interact according to a reduced Hamiltonian of the form 

H = lkr n = E V & ~ °i) > W 

<ij> 

the sum running over nearest-neighbour pairs of sites < ij >. Imposing V(ff) = 
V(—0), a particular system is characterised by the three numbers V r = V(irr/3) — 
V(0), r = 1,2,3. A duality transformation maps the space of such models to itself 
0-01 : with the couplings parametrised by variables x r = exp(— V r ), a sum over the 9i 
is equivalent to one over dual variables 6$, with the triplet of couplings {x r } replaced 
by the dual set {x r }: 

X\ = (1 + X\ — X 2 -X 3 )/A 
X 2 = (1 — X\ — X 2 + X 3 )/A 

x 3 = (1 - 2xi + 2x 2 - x 3 )/A (2) 
where A = 1 + 2x% + 2x2 + #3- The transformation leaves invariant points on the line 

(xi, X2, x 3 ) = (t, [3— at, a— 2/3i) (3) 

where a = 3— V6 and (3 = v6— 2. This self-dual line intersects the line of 6-state 
Potts models at the point x\ = X2 = X3 = (v6— 1)/5, denoted P below. There is also 
a distinguished surface on which the models decouple into independent three-state 
Potts and Ising models H. This is revealed on writing each spin Si as a product of 
a Z2 and a Z 3 valued variable: 

Si = e"* = E i = l,e ±2w< / 3 J a = ±1 (4) 

The pairwise interaction energy can then be written as 

V(0) = Ji[l- cos(0)] + J 2 [l- cos(2#)] + J 3 [l- cos(3#)] 

= Jill-i^+.S- 1 )] + Jsfl-^S+S- 1 )] + J 3 [l-±(a+0] 

= 3Ji[1-^ iCT / EiEj ] + l(J 2 -J x )[l-8^] + (2J 3 -J 1 )[l-6 ai(r .} (5) 

where 6 = 0i—6j, S = S{/Sj, £ = £j/£j, a = (Ti/cTj, and 

j x = 1 In _a_ j 2 = 1 l n JEa_ j 3 = 1 i n 4L . ( 6 ) 

The spins £j and cr, decouple on the surface Ji = 0, xi = x 2 x 3 , on which I J2 is the 
3-state Potts coupling, and 2 J 3 the Ising coupling. 

The more general phase structure of these models has been studied by various 
authors over the years (see for example [0-0] ) . The question is an interesting exercise 
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in its own right, and is also of wider relevance - to, for example, the effect of hexag- 
onal symmetry-breaking on the isotropic planar model |§,|6|, the behaviour of the 
cubic model 0,0], and the spectra of Heisenberg antiferromagnetic spin chains ||. 
An initial phase diagram was suggested by Domany and Riedel in 0, with the 
three-dimensional thermodynamic phase space partitioned into four domains, one 
disordered and the others exhibiting Zg, Z3 and ordering. The existence of an 
additional massless phase was then demonstrated, first j|] along the Villain [ 10 1 line, 
and then || throughout a whole three-dimensional region. This phase was incor- 
porated into the Domany-Riedel diagram by Alcaraz and Koberle in Q, its end 
point on the self-dual line later being identified with a particular integrable point C 
on the phase diagram (further support for such an identification can be found 

EM)- 
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However there are reasons to believe that the story is not yet complete. We give 
two examples. 

(i) In the diagrams of [Q] and Q , the point P touches surfaces of transition from the 
disordered phase into regions of Z2 and Z3 order. Such transitions are expected to 
be of second order, and so the correlation length on these surfaces should be infinite. 
This contradicts the known behaviour of the 6-state Potts model at P, where the 
correlation length remains finite [14], albeit large [15|. 

(ii) On the decoupling surface J\ = 0, there is a line of Ising transitions, J3 = 
^ ln[l+\/2], and a line of 3-state Potts transitions, J2 = | ln[l+\/3]- These lines 
cross at a renormalisation group fixed point D, the product of a critical Ising model 
and a critical three-state Potts model. Three Zg-invariant operators at this point 
are the Ising and three-state Potts energy densities e and £, and their product e£. 
Since their scaling dimensions are 1,4/5 and 9/5 respectively, all are relevant and the 
fixed point is triply unstable. This corrects the approximate (Migdal) renormalisation 
group result used in , which gave D as being doubly-unstable. 

We can take this second point a little further, using continuum field theory ar- 
guments which are valid in the scaling region around the point D. The fixed point 
itself is described by a product of c = 1/2 and c = 4/5 conformal field theories, and 
nearby points by perturbations of this product c = 13/10 conformal field theory by 
combinations of the continuum operators e, £ and e£. The first two are anti-self-dual 
under @, while the third is self-dual and moves the model away from D along the 
line (||). Minimal models coupled by local operators have received a fair amount 
of attention (see, for example, ]|1|) but this particular instance does not seem to 
have been studied in any detail. However, Zamolodchikov's counting argument [17] 
can be used to show that the e£ perturbation preserves at the very least conserved 
charges of spins ±3 and ±5. The perturbed (continuum) theory should therefore be 
integrable [J^] and we can hope to obtain information about the scaling region of the 
self-dual line near to D via the thermodynamic Bethe ansatz (TBA) technique [19]. 
This method expresses the finite- volume 'effective central charge' c(r), r = mR, of 
a model with bulk length scale 1/m confined to a circle of circumference R in terms 
of the solutions e a (9) to a set of coupled integral equations. A candidate system for 
this case has already been identified: in [20] it was observed that the following TBA 
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system for the functions e\ . . . £7 



e a {9) = 5 al r coshO - — Y,l [ ab ] ^*H^+e- eb W) , 

b=i 

o roc 

c(r) = — d9r cosh 9 ln(l+ e - £l(e) ) 



(7) 



predicts an ultraviolet central charge of 13/10, and a small-r behaviour of c(r) 
compatible with an expansion in even powers of a coupling A to an operator of 
dimension 9/5. (Here * denotes the convolution, f*g(9) = dB r f(6—9')g(6'), 

4>(9) = l/cosh#, and is the incidence matrix of the Ej Dynkin diagram, with 
1 labeling the node at the end of the middle- length arm.) This is the behaviour 
expected of the perturbation of the product Ising and three-state Potts conformal 
field theory by its e£ operator. (On dimensional grounds, m will be related to A as 
m oc A 5 , but we do not attempt to find the constant of proportionality here.) Assum- 
ing that the TBA is correct, it can be used to extract a non-trivial prediction about 
the vacuum structure of the perturbed model, using arguments described in section 
4.3 of jHJ. The system (^) implies the following asymptotic for the ground-state 
energy E(m, R): 



TT 



E(m,R) = c(mR) 



m 



TT 



K\ (mR) 



(8) 



where K\ is the modified Bessel function of order one and the prefactor v6 can 
be interpreted as the Perron-Frobenius eigenvalue of the incidence matrix for single 
kinks interpolating degenerate vacua of the perturbed theory fi~3|| . These vacua must 
support a representation of the global 7L§ symmetry, and we find that the simplest 
incidence matrices compatible with both this fact and the eigenvalue implied by @ 
are the following: 
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(9) 



The first is consistent with the coexistence of disordered and ordered phases, the 
second with coexistence of Z2 and Z3 ordered phases. This leads us to the first part 
of our proposal for the resolution of points (i) and (ii) above: we suggest that the 
entire segment of the self-dual line from D to C, which includes the 6-state Potts 
transition point P, lies on a surface of first-order transitions separating disordered 
and fully-ordered regions, and in the scaling region near to D is described by an 
integrable scattering theory with kink incidence matrix I a . (In the massive scaling 
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region near to C, the vacuum structure comes as the iV=6 case of the results of [13], 
and is also given by I a .) It is natural to suppose that the matrix I b describes the 
vacuum structure on the opposite side of the point D, and support for this idea 
comes from the following argument. 

Consider the behaviour of the model near to the decoupling surface J\ = 0, 
taking either J 2 S> 1 or J3 3> 1. In the first case, the effect will be to freeze out the 
three-state Potts spins "Si, and all terms can be approximated by 1 in the final 

line of (|5|). The result is an effective interaction for the remaining unfrozen spins, 
equal to 2(J3+Ji)[l— 5 ai(Tj \. Thus the only effect of a non-zero J% in this 'E-frozen' 
region is to replace J3 with J3+J1. The line of Ising transitions on the decoupling 
surface was at J 3 = J 3 C = \ ln[l+V2] ; we now conclude that a small non-zero J\ in 
the region J2 3> 1 simply shifts this line to J3 = J|— J\. A similar argument shows 
that the line of three-state Potts transitions, situated at J 2 = J 2 = | ln[l+v3] on 
the decoupling surface, is shifted to J 2 = J 2 — Ji i n the region J3 S> 1 where the a 
spins are frozen. Finally, duality can be used to see that in the opposite regimes 
J2 -C 1 and J3 -C 1, the critical lines move in the opposite senses when J\ shifts 
away from zero. This line of argument has nothing to say directly about the central 
region where the Ising and three-state Potts spins are both near their critical points, 
but the simplest hypothesis, also consistent with the continuum results of the last 
paragraph, is to continue the lines in from the asymptotic regions as in Fig. 1. 




Figure 1: Schematic pictures of slices through the phase diagram just below, 
and just above, the decoupling surface. 

Much of our numerical work has been devoted to the confirmation of this picture. 
Before we describe this, we pause to introduce a third set of coordinates on the 
thermodynamic phase space. First, notice that the parameter K defined by 

K(x 1 ,x 2 ,x 3 ) = (x 1 -x 2 x 3 )/A (10) 

is mapped into itself by duality: K = K . Therefore surfaces of constant K in phase 
space are mapped into themselves by duality; K = is the decoupling surface. Points 
with K > lie 'below' this surface, on the same side as the 6-state Potts point P. 
It is easy to check that the two parameters 

1 + 2X2 , I+X3 
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have the simple duality transformation rule 



m = ~V2 m = -3/3 (12) 

It is convenient to use (K, 2/2,2/3) as coordinates in phase space: on any fixed- K 
surface the point (K, 0, 0) is the self-dual point and duality is just the reflection in 
this point. The origin (0, 0, 0) is the fixed point D on the decoupling surface, and 
the 6-state Potts transition point P is at (K P ,0,Q) with K P = (^/6-l) 2 /25 « 0.084. 
On the decoupling surface, y% and 2/3 are functions of the Ising and Potts couplings 
respectively. 

A cluster algorithm appears to be the ideal choice to simulate the relevant regions 
of the phase diagram, since it is reasonable to expect large correlation lengths even 
when the transitions are only first-order. Since cluster algorithms are especially 
simple to implement for Z2 and Z3 models, we exploited the possibility of writing 
the Hamiltonian as in (|5|): the algorithm we used performs alternate cluster updates 
of the o"j and Sj variables, where those which are not being updated provide effective, 
site-dependent couplings for the others. For example, suppose we are updating the 
Z2 variables: the effective coupling to be used on the link < ij > is then J3 + 
J\ (E+S _1 ) /2 with £ = Sj/Ej, while when updating the Z3 variables the effective 



coupling is J 2 + J\o with a = a-i/aj. Similar algorithms were introduced in [21] for 
the Ashkin Teller model, defined as two coupled Ising models. 

The fact that the effective couplings are site-dependent does not pose a problem 
for the cluster updates. However the effective Z2 and Z3 couplings can become nega- 
tive, i.e. antiferromagnetic. This will happen for the effective Z2 coupling whenever 
^3 + ^1 < or J3 — Ji/2 < 0, and for the Z3 coupling whenever J2 + J\ < or 
J2 — J\ < 0. Negative couplings on some links can in turn lead to frustrations, which 
in principle can make the cluster algorithm highly ineffective. However it is easy to 
convince oneself that the Z3 model is never actually frustrated: for every configu- 
ration of the Z2 spins there exists a configuration of the Z3 ones such that all the 
links are satisfied. Therefore the possibility of frustrations exists only for J\ < — J3 
or Ji > 2J3. These relations are never satisfied in the regions we considered. 

We explored several fixed- K surfaces and mapped the various phases and their 
boundaries using the Binder cumulants method In our case, we define two 

cumulants, one for each order parameter: 

_3 1 (mi) (m|) 
Q °~2~2^ ^- 2 "^|)2 ( 13 ) 

where m a and are the Z2 and Z3 magnetizations per site. For each surface we 
measured the Z2 and Z3 Binder cumulants on a grid of points around the self-dual 
point (K, 0, 0) for two different lattice sizes L\ x L\ and L2 XL2. For every fixed value 
of y2 we estimated the value of j/3 where each Binder cumulant remains constant when 
the lattice size is increased, i.e. the transition point. This allows us to determine a 
set of points belonging to the transition lines for the Z2 and Z3 variables. These are 
the points plotted in Figs. 2,3,4 for the K = 0.06, 0.04, —0.08 surfaces respectively. 
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The consistency of the procedure was checked by repeating it with the roles of and 
ys inverted. Finite-size effects are signaled by violations of the duality symmetry, 
that is, in our coordinates, symmetry under reflection in the origin. Small violations 
are visible in our figures; we checked that they become smaller when the lattice sizes 
are increased. These effects turned out to be much larger in the K < region. In 
fact we observed that to control them satisfactorily in the K > region lattice sizes 
as small as L\ = 15, L 2 = 30 sufficed, while for K < we had to use L\ = 45, 
L 2 = 60. 

The numerical results confirm the qualitative picture shown in Fig. 1, and our 
claim that the decoupling point D on the self-dual line marks a change from the 
coexistence of disorder and order to the coexistence of Z2 and Z3 order. The 
reduced length of the first-order segment at K = 0.04 is consistent with scaling 
predictions, though we are not close enough to the point D to see a complete collapse 
of data. In conclusion, we have proposed a significant modification to the previously- 
accepted phase diagram of the spin models, and this has been supported by a 
detailed numerical study. 



0.06 



0.04- 



0.02- 




-0.02 



-0.04 



-0.06 



-0.04 -0.03 -0.02 -0.01 0.01 0.02 0.03 0.04 

Y2 

Figure 2: Phase diagram on the K = 0.06 surface. Phases are labelled accord- 
ing to the nature of their ordering. Binder cumulants for lattice sizes L\ = 15 
and L 2 = 30 were used. 
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Figure 3: Same as Fig. 2 for K = 0.04 
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Figure 4: Same as Fig. 2 for K = —0.08. The lattice sizes used here were 
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